Dissipative Preparation of Spin Squeezed Atomic Ensembles in a Steady State 
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We present and analyze a new approach for the generation of atomic spin squeezed states. Our 
method involves the collective coupling of an atomic ensemble to a decaying mode of an open optical 
cavity. We demonstrate the existence of a collective atomic dark-state, decoupled from the radiation 
field. By explicitly constructing this state we find that it can feature spin squeezing bounded only 
by the Heisenberg limit. We show that such dark states can be deterministically prepared via 
dissipative means, thus turning dissipation into a resource for entanglement. The scaling of the 
phase sensitivity taking realistic imperfections into account is discussed. 
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The realization of spin squeezed states [1 of atomic 
ensembles is an important subject in quantum science. 
Such states play a central role in studies of many-body 
entanglement [2]-[4]. In addition, they may lead to the 
practical improvements of state-of-the-art atomic clocks 
and frequency standards ^5ttZj- In a spin-squeezed-state 
the atoms are entangled in such a way that the fluctua- 
tions of their total spin are smaller than the sum of the 
fluctuations of the individual atoms. As a rule, such en- 
tangled states are extremely fragile and spin squeezing is 
destroyed due to dissipation or decoherence. In this Let- 
ter we propose a new approach for the realization of spin 
squeezing using recently developed ideas on quantum- 
bath engineering [8UT2] and show how to realize squeez- 
ing in the steady state of a dissipative atom-cavity sys- 
tem. As a result, our approach can be considerably more 
robust against noise as compared to existing preparation 
methods [13-19 that produce short-lived spin squeezed 
states, limited by decoherence processes. Moreover, using 
properly selected atomic transitions allows to continously 
pump the atoms into the desired states and consequently 
avoid population losses due to scattering into states not 
not participating in the squeezing. A similar approach 
to continously entangle two remote atomic ensemble us- 
ing free space scattering has been experimentally demon- 
strated in [20j. 

The central idea of our work can be understood by 
considering an ensemble of N atoms interacting with a 
single radiation mode of an open optical cavity (Fig. la) 
and externally driven by a pair of coherent laser fields. 
The cavity mode and laser fields are tuned to excite a 
pair of two-photon Raman transitions, each involving one 
laser field with Rabi frequency and a single cavity 
photon a {a^) (Fig. lb). Assuming that the coupling of 
atoms to the cavity mode is uniform, the unitary time- 
evolution of such cavity-atom system is described by the 
effective Hamiltonian 

i/int = |-at ((^+5+ + ^-S-) + H.c, (1) 
where is the spin raising (lowering) operator 




FIG. 1. (a) Schematic experimental setup of atoms with a 
spontaneous emission rate 7 in a cavity with line-width Hi 
supporting a single cavity mode a. (b) Effective linkage pat- 
tern consisting of two degenerate ground states ( | ± ) ) en- 
coding the effective atomic spin and two excited state \e±) . 
The system is driven via classical control fields Q± with large 
detuning A and a single global cavity mode a coupling homo- 
geneously to all atoms. 



of the collective spin, g is the single atom-field cou- 
pling strengthand A is the single-photon detuning from 
the excited state. The coupling gQ±/A is therefore an 
effective Rabi frequency corresponding to off-resonant 
Raman scattering (cf. Fig. [T]). Direct examination 
of the Hamiltonian ([T]) shows that it has a dark-state 
I ^ ) = I V^spin ) I Ocav ) , where | 

cav ) is the cavity vac- 
uum and I V^spin ) satisfies 

{n+s+ + n_s-) li^spin) =0. (2) 

The interesting features of this state are (i) that it is a 
pure state not containing any photons and consequently 
is not affected by the cavity decay, i.e. | V^spin ) is a dark- 
state of the cavity decay, and (ii) it describes a highly 
correlated atomic state containing a high degree of spin 
squeezing. The latter can be seen by computing the ratio 
between the spin fluctuations in the x and y direction: 

^-\-)^ / +r^+)^. For any 7^ this ratio is smaller 
than one, indicating that the dark state is spin squeezed. 

To quantify the degree of useful spin squeezing, we 
use the phase sensitivity introduced by Wineland et al. 
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FIG. 2. Phase sensitivity of the dark- state | V^spin ) : exact nu- 
merical solution (solid curve) and mean-field approximation 
(dashed curve) for = 10 (left) and N = 100 (right) atoms. 
The optimal squeezing is obtained for r2+ Q- and ap- 
proaches the Heisenberg hmit Scj) = l/y^N{N/2 + 1) ^ 1/N. 
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We first focus on symmetric states with the total spin 
S = N/2, i.e. an initially spin-polarized ensemble, and 
expand the dark-state into eigenstates of Sz as | V^spin ) = 
Em Cm I ^ = N/2, Sz = -N/2 + m ) . Substituting this 
into eq. Q results in a recursion relation between Cm 
and Cm-\-2^ in close analogy with squeezed light [22 
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Here (:) are the binomial coefficients and cq is deter- 
mined by the normalization condition Em = 1- 
The phase sensitivity S(j) of | V^spin ) is shown by the blue 
solid curve in Fig. |2]as function of the ratio Q-^/Q- be- 
tween the control fields. For Q_^/Q_ = 0, the dark-state 
is fully polarized along the z-axis, giving a phase sensi- 
tivity of l/\fN corresponding to the Standard Quantum 
Limit (SQL). In the opposite limit of Vtj^ the 
phase sensitivity approaches the Heisenberg limit [151 ES] 
S(l) = l/^N{N/2 + 1), indicating that the dark state 
corresponds to an almost maximally squeezed atomic 
state. 

We now turn to the preparation of the dark state ([T]) 
and study the corresponding quantum dynamics, focus- 
ing first on the effect of the cavity decay. To this end, 
consider an ensemble of N four-state atoms, consisting of 
two (meta-) stable states | + ) , | — ) serving as our effec- 
tive spin-states and two excited states \e±) . The atoms 
are coupled to a single mode of an open cavity (Figjl]) 
with resonance frequency cj^, volume V and photon es- 
cape rate n driving the transitions | ± ) \^±) and 
strong control fields Vt± with frequency uJc driving the 
transitions | ^ ) I e± ) , respectively. The atom-cavity 
field coupling strength is given by ^ = p / h^/hIJZj2e^ , 
where p is the dipole matrix element of the transition. 
We assume the control fields and the cavity mode to 
be far detuned from the resonant transition frequency 
cjres = ^± — ^e± between ground state and excited state 



manifold. This allows us to adiabatically eliminate the 
excited states and we obtain the effective Hamiltonian 

§ with s+ = Ef = Ef 1 + ),, (-1. = 

and |A| = \uJres — ^o,c| El- Assuming hi ^ |A|, we 
can neglect the occupation of the cavity mode and using 
standard techniques [25] obtain the master equation in 
Lindblad form after tracing over the cavity modes 
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{XtX,p}^-2XpXt 



(5) 



Here {•,•}+ denotes the ant i- commutator, jcav = 
g^n^/A^K with = nl^nl, I = sin OS^ + cos6>5'_ 
is the Lindblad-operator accompanying the emission of a 
cavity photon, and we defined tan^ = 

To compute the phase sensitivity S(j), we use (§ to 
determine the time evolution of the expectation values of 
the spin operators Sz and S^ 



d{Sz) _ 

dt 
d{Sl) _ 

dt 



-7cav [COS\0){S+S-) 

-7cav(sin(6>) - cos((9)) 



^\e){s-s^)] (6) 

(7) 

{{sm{e)S- + cos(^)5+) {S:,Sz + SzS:,) + H.c.) 

Note that the dark state | V^spin ) , defined by p|, is in- 
deed a steady state of ([6][7|. Equations (|6]|7| involve 
higher moments of the spin, whose time evolution should 
be computed independently. To avoid this complication 
and obtain an analytic solution, we linearize the equa- 
tions of motion around {Sz) ^ —N/2. Defining the small 
fluctuations as 5Sz = 5^ + N/2 and using the approxi- 
mation {Si + Si) = (S^) - {Si) ^ N{5Sz) + N/2 we find 
after some algebra 



dt 
dt 



-TcaviV [cosi2e){SS,) - sin^e)] , (8) 

1" 



cos(26')(S'^) + (l-sin(26i))- 



.(9) 



The solutions to these equations decay exponentially in 
time with an effective rate 7eff = 7cav cos{20)N. From 
eqs. ([8][9| we obtain the phase sensitivity in the steady 
state as 



7V(1 -sin(2(9))cos(2i9) 
(7Vcos(26>) -2sin2((9))^ 



(10) 



Looking at eq. ( |1Q[ ) it is obvious that 1]+ ^ ft-, or 
equivalently 6 ^ 7r/4, is necessary to obtain a high de- 
gree of spin squeezing. However, in this limit the effec- 
tive dark-state pumping rate 7eff tends to zero making 
the spin squeezing process extremely slow. This renders 
competing processes, as e.g. spontaneous Raman scatter- 
ing of individual atoms, very important as we will see be- 
low. Moreover, as can be seen in Figj2j this approximate 
solution (dashed red curve) agrees well with the exact 
result (solid blue curve). As expected, the two curves 
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deviate significantly only for 6 = tan ^ 
7r/4 — 0{1/N)^ where the linearization fails because 
{6S,) = 0{N). 

Up to this point, we have restricted ourselves to the 
manifold of total spin S = N/2. This approach is valid 
if the system is initially prepared in the maximally po- 
larized state I S = 5'^ = N/2 ) because Hamiltonian 
([l]) commutes with the total spin S^. Consequently, if the 
system is initially prepared in a state with a different to- 
tal spin, the final steady state will be different and will 
in general contain a smaller amount of spin squeezing. 
This is shown in Fig. [sja), where we compare the lin- 
earized solution with the exact numerical solution of the 
Lindblad eq. ([5| for = 10 spins with different initial 
conditions. The linearization works well only when the 
spins are initially polarized, as would be expected. 

So far our analysis completely disregarded the effect of 
spontaneous Raman scattering of the individual atoms 
into free space. Being a single atom process, it breaks 
the conservation of the total spin and thus competes with 
the squeezing process. To investigate the effect of sponta- 
neous Raman scattering we introduce a spontaneous de- 
cay rate 7 <C |A| to the states \e±) . The corresponding 
master equation is easily obtained using standard tech- 
niques and reads 



d 



(11) 



N p 

Tspont ^ ^ ^ ^ |^>^j,Q!^>^j,Q!7 Pj* '^^j^aP^j^a^ 



The first term on the RHS of (11) is equivalent to eq. 
([5|, whereas the second term describes the spontaneous 
Raman scattering of atom j into free space with rates 
^Vfi / /S? . The Lindblad-operators are given by 



Tspont 

sin 6>cr 



(cos ^(7+ + singer and C 



The numerical solution of (11) is shown in 
for different initial conditions and = 10 
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Fig. ^ 

spins. The spontaneous Raman scattering couples differ- 
ent total spin-S* manifolds and, as expected, lifts the de- 
generacy between them. Consequently it leads to a mixed 
but nevertheless unique steady state, characterized by a 
unique density matrix, independent of the initial state, 
and hence presents rather a feature than a detrimental 
problem. 

Still, Raman scattering leads to a reduction of the 
achievable phase sensitivity. To quantify this effect we 
obtain the contribution of the single-spin decay to the 
time evolution in the linearized regime by using the sec- 
ond row of ( pTj ) to derive 



dt 
d{Sl) 
dt 



Tspont 



= -27, 
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FIG. 3. Time evolution of the phase sensitivity for N — 10 
spins and = 0.2. Numerical solution of the mas- 

ter equation, with fully polarized (solid curves) or randomly 
generated (dashed-dotted curves) initial state, and mean- field 
approximation (dashed curves), (a) The atoms are coupled 
only to the cavity, i.e. 7 = 0. (b) Finite scattering rate 
7 7^ with single- atom cooperativity x = 1- Observe that 
the steady state does not depend on the initial state. 



to the cavity, eqs. ([8j[9|), we find that in the steady state 
the competition between the two processes leads to: 



^cos(20)+7 ' 
1 (1 _ sm{29)) + 27 (1 - 1 sin(2^)) 
4 ^ cos(2^) + 27 (1 - i sin(2^)) 



(14) 



(.15) 



Introducing the single-atom cooperativity x = 
Tcav/Tspont = 9^ / ^1 ^ the optimal value of the phase 
sensitivity (50 can be obtained by numerically minimizing 
the corresponding expression with respect to the mixing 
angle d while keeping x fixed. The resulting curve is 
shown in Fig. [|for A/- = 10^. To gain insight into this 
result, we will now separately investigate the limits of 
small and large single-atom cooperativity. 

For small cooperativity, i.e. ^ I7 optimal spin 



squeezing is obtained for d <C 7r/4. Expanding eqs. (14 
[l5| around ^ we obtain a phase sensitivity 5(jr~~^ 
(1 — A^x^ + 46>^) with the optimum value obtained for 
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Adding to this the contribution of the collective coupling 



It should be noted that due to the collective enhancement 
by a factor of A/", one gains a quadratic improvement in 
the phase sensitivity even for small cooperativities. 

For a large cooperativity X 1? optimal squeezing 
is obtained for 17+ ^ and one may naively expect 
the linearization to fail. However, to render the problem 
tractable we nevertheless perform the large- expansion 
and explicitly check the validity of the obtained results 
below. To obtain the scaling of the phase sensitivity we 
expand eqs. (14][l5) around = 7r/4 — e. For small 
e ~ 1/ a/TV <C 1 and neglecting terms of order 1/N we find 
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FIG. 4. Optimal phase sensitivity in the steady state, in the 
presence of single spin decay for N = 10^. The solid curve 
is given by the numerical minimization of e qs.(|14||15 |). The 
dashed-dotted and dashed curves are eqs. ( |16|17[ ), 
tively. 
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2 _ X i+2A^x€ which has a minimum at e = ^/l/2NX' 
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This leads to a sensitivity 
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To check the validity of the linearization, we compute 
{n)/N for the optimal value of and find {n)/N = 
^j^t^^ ^ 1/a/]V <C 1, validating its use in the limit of a 



large ensemble. Eq. (17) shows that the achievable phase 
sensitivity scales as S(j) ^ N~^^^ and thus offers a signif- 
icant improvement over the SQL. We emphasize that in 
our case spin squeezing is deterministically generated in 
a steady state and can thus potentially be more robust 
against external perturbations as compared to earlier ap- 
proaches p!5l [26]. 

A possible experimental realization could be set up 
in ^^Rb using the clock states |F = l,mi? = 0) and 
\F = 2^mF = 0) in the 551/2 ground-state manifold. 
Using, e.g.. circularly cr+-polarized cavity and control 
fields, couples these states to the | F = 1, = +1 ) and 
I F = 2, = +1 ) of the 5Pi/2-manifold. Using param- 
eters of current experiments p!7|, we have N ^ 10^ and 
X ^ 0.1, predicting a phase sensitivity S(j)/5(j)sQL ^ 0.07, 
i.e. an improvement of more than one order of magnitude 
with respect to the SQL. 

Apart from the spontaneous Raman scattering, there 
are two more effects that have been neglected. First, 
spontaneous Raman scattering in a multi-level atom, 
such as the suggested ^^Rb, can also lead to states out- 
side of the effective four-state system. However, as shown 
above, the squeezed steady state is (mixed but) unique 
and is achieved from any intial state. This allows us to 
use additional optical pumping fields that repump the 
atoms into the correct configuration. In the case of ^^Rb 
these would be linearly polarized fields resonantly driv- 
ing the F = 1 ^ = 1, F = 2 ^ = 2 and Amp = 
transition, that leave the clock states unaffected as these 
transitions are forbidden. The repump-processes can eas- 
ily be incorporated in our formalism as additional single 
atom scattering channels resulting in an effectively larger 



Raman scattering rate. Consequently only the coopera- 
tivity X is reduced and our scaling results are not affected. 

Second, we assumed a spatially homogeneous coupling 
of the atoms to the cavity mode, which can only be ap- 
proximately achieved in a ring cavity. For a standing- 
wave cavity, the coupling-constants depend on the po- 
sition of the atoms, and consequently change in both 
absolute value and phase. However, variations in the 
magnitude of gj cause only a small reduction of the ef- 
fective spin cooperativity • This can be understood by 
observing that amplitude variations mainly preserve the 
symmetry of the collective state. On the contrary, fluctu- 
ations in the phase of the single atom coupling constants 
gj are detrimental as they effectively lead to dephasing 
of the spin squeezed state and lead out of the symmetric 
subspace. To avoid this effect, it is necessary to fix the 
atoms in space, which can be achieved using standard 
cooling and trapping techniques. 

To summarize, we have presented a method for the de- 
terministic generation of spin squeezed states in a driven 
ensemble of effective two-state atoms in a strongly dissi- 
pative cavity. The generating process can be understood 
as optical pumping into a non-equilibrium steady state 
that at the same time is a dark-state with respect to the 
cavity decay. Introducing spontaneous Raman scattering 
we showed that the squeezed steady state is unique and 
does not depend on the initial state of the system. We 
discussed the effect of the single atom spontaneous Ra- 
man scattering on the achievable phase-sensitivity and 
found that it scales favorably with the single-atom coop- 
erativity, indicating that the present method can be of 
direct importance for, e.g., optical atomic clocks. 
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